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Abstract. Recently we have evaluated the matrix elements (Or p ), where O — {1, 0, zan/J} are the 
standard Dirac matrix operators and the angular brackets denote the quantum-mechanical average 
| for the relativistic Coulomb problem, in terms of generalized hypergeometric functions 3F2 (1) for 

all suitable powers and established two sets of Pasternack-type matrix identities for these integrals. 
\ The corresponding Kramers-Pasternack type three-term vector recurrence relations are derived here. 
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1. Introduction 



Recent experimental and theoretical progress has renewed interest in quantum electrodynamics of 
atomic hydrogenlike systems (see, for example, [22], [23], [27], [28], [36], [53], and [55] and references 
therein). In the last decade, the two-time Green's function method of deriving formal expressions 
for the energy shift of a bound-state level of high-Z few-electron systems was developed [53J and 
numerical calculations of QED effects in heavy ions were performed with excellent agreement to 
current experimental data [22] , |23j . [55] . These advances motivate, among other technical things, 
investigation of the expectation values of the Dirac matrix operators between the bound-state 
relativistic Coulomb wave functions. Special cases appear in calculations of the magnetic dipole 
. hyperfine splitting, the electric quadrupole hyperfme splitting, the anomalous Zeeman effect, and 
the relativistic recoil corrections in hydrogenlike ions (see, for example, [I], [52], [51], [57] and 
references therein). 

as : 

In the previous paper [57], the matrix elements (Or p ), where O = {1, (3, ian/3} are the standard 
Dirac matrix operators and the angular brackets denote the quantum-mechanical average for the rel- 
^ ativistic Coulomb problem, have been evaluated as sums of three special generalized hypergeometric 
functions 3F2 (1) (or Chebyshev polynomials of a discrete variable) for all suitable powers. As a 
result, two sets of Pasternack-type matrix symmetry relations for these integrals, when p — > —p — 1 
and p — > —p — 3, have been derived. We concentrate on what are essentially radial integrals since, 
for problems involving spherical symmetry, one can reduce all expectation values to radial integrals 
by use of the properties of angular momentum. 

Our next goal is to establish relativistic analogs of the Kramers-Pasternack recurrence relation 
[3~T] . |44j . and [15] (some progress in this direction had been made in [2]). Here several three-term 
vector recurrence relations are obtained, which follow in a natural way from the well-known theory 
of classical orthogonal polynomials of a discrete variable [39J and basic matrix algebra. This paper 
is organized as follows. We review the nonrelativistic case in the next section and then present 
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relativistic extensions, that seem to be new, in sections 3 and 4. The last section contains a new 
interpretation of the known two-term recurrence relations for these relativistic expectation values 
[50] . |51j . The new recurrence relations are summarized in appendix A in matrix form for the benefit 
of the reader. 

The author believes that these mathematical results are not only natural and elegant but also 
will be useful in the current theory of hydrogenlike heavy ions and other exotic relativistic Coulomb 
systems. At the same time, the mathematical structure behind these expectation values remains 
not fully understood [2] . For example, to the best of my knowledge, numerous recurrence relations 
for the radial integrals, obtained with the help of a hypervirial theorem (see PQ, [50], [51], [51] 
and references therein), do not follow in an obvious way from the advanced theory of generalized 
hypergeometric functions [5J. After more than 80 years of a thorough investigation, the relativistic 
Coulomb problem keeps generating some mathematical challenges. 



2. Expectation Values (r p ) for Nonrelativistic Coulomb Problem 

Evaluation of the matrix elements (r p ) between nonrelativistic bound-state hydrogenlike wave 
functions has a long history in quantum mechanics. An incomplete list of references includes 

0, 0, 0, is], 0, ma, pi us, pi, pi, ubi pi, una, us], [33, isg, 01, isa, [33, us, 

H, @S], HZ], US], ED], [2S], [SB], [2S], [52], and [BJ. Different methods were used in order to 
derive these matrix elements. For example, in Ref. [58] (see also [H], [15]) the mean values for 
states of definite energy have been obtained in terms of the Chebyshev polynomials of a discrete 
variable t m (x,N) = hm °' ) (x,N) originally introduced in Refs. [60] and [6i] . The so-called Hahn 
polynomials, introduced also by P. L. Chebyshev [62J and given by 

were rediscovered and generalized in the late 1940s by W. Hahn (we use the standard definition of 
the generalized hypergeometric series throughout the paper [6J, [2~T]). 

The final results have the following closed forms 

^-'H^Gfr^-'- 1 -- 2 '- 1 )' < 2 - 2 > 



when k = 0,1, 2, ... and 



1 / 17. x //+2 



t fc (n-Z-l,-2Z-l), (2.3) 



^ r k+2 J 2n \nao 

when = 0, 1, 2Z. Here clq = h 2 /me 2 is the Bohr radius (more details can be found in Refs. 
and [58]). 



The ease of handling of these matrix elements for the discrete levels is greatly increased if use is 
made of the known properties of these classical polynomials of Chebyshev [2T] , [26] . [29], [39], [4"0] 
and [60], [61], [62]. The direct consequences of fl2.2l) - fl2.3l) are an inversion relation: 

2k+l (0J _ UM 

(r*- 1 ) (2.4) 



1 \ (2Z\ 2k+1 (21 - k)\ , A _ h 



r k+2 I \na J (21 + k + 1)\ 
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with < k < 21 and the three-term recurrence relation: 



^ = 2n (2fc + l) /na \ , k _^ _ k ((21 + l) 2 - k 2 ) /nao\2 2 ^ 



(2.5) 



+ 

with the initial conditions 

-) = A> (!> = !. (2-6) 
which is convenient for evaluation of the mean values (r fc )> for k > 1. 

In our approach, the recurrence relation (12.51) is a special case of the three-term recurrence 
relation for the Hahn polynomials (13.81) . It was originally found by Kramers and Pasternack in 
the late 1930s [31], [33], and [15]. The inversion symmetry (12.41) . which is also due to Pasternack, 
has been rediscovered many years later [12], [37] (see also [23] and [25] for historical comments). 
Generalizations of (I2.4l) - (l2.5p for off-diagonal matrix elements are discussed in Refs. [TTj . [T7] . 
ED], [25J, [3E], 01] , [32], [3E], [31] , [50], [51], [5S], and [53]. The properties of the hydrogenlike 
radial matrix elements are considered from a group-theoretical viewpoint in Refs. [3], [13], and 
[32]- Extensions to the relativistic case are given in [TJ, [2], [IB], [TjJ], [33], [57], and [5B] (see also 
references therein and the following sections of this paper). 



In a retrospect, Pasternack's papers [33], [32] had paved the road to the discovery of the continuous 
Hahn polynomials in the mid 1980s (see [3], [5J, [30], [58] and references therein). 

3. Three-Term Recurrence Relations for Relativistic Matrix Elements 

In the present paper we establish relativistic analogs of the Kramers-Pasternack recurrence rela- 
tions (12.51) for the following set of integrals of the relativistic radial functions: 

POO 

A p = / r p+2 (F 2 (r) + G 2 (r)) dr, (3.1) 
Jo 

POO 

B p = / r p+2 (F 2 (r) - G 2 (r)) dr, (3.2) 
Jo 

POO 

C p = r p+2 F(r)G(r) dr. (3.3) 



Jo 

With the notations from Refs. [57] and [58], one can evaluate these integrals in terms of the Cheby- 
shev and Hahn polynomials of a discrete variable and present the answer, say, when p > 0, in the 
following closed form: 

Afiu 2 {2a(3) p A p (3.4) 
= clk (en + v) h p ^ (n — 1, —1 — 2u) 

_2^±J^ (s 2 k 2 - v 2 ) h£> 1] (n - 1, -1 - 2u) 
p+1 K j P \ ) 

+an (sk — v) h p \^ (n, 1 — 2u) , 

Afxu(2ap) p B p (3.5) 
= a (sk + v) h p \^ (n — 1, —1 — 2v) 
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-a (en — v) h p °l^ (n, 1 — 2v) , 



^v 2 (2aP) p C p 



(3.6) 



a/j, (en + v) h p °+^ (n — 1, —1 — 2 
_ 2 ^±^ ( £ V - u 2 ) hp 1 ' V (ji _ i, 



p + 



2v) 



+a/i {en — v) h p °+^ (n, 1 — 2v) . 



Here, 



K = ±{j + 1/2) , v = v 7 ^ 2 - M 2 , 
/x = aZ = Ze 2 /hc, a = y/l — e 2 , 
e = E/mc 2 , (3 = mc/h 
with the total angular momentum j = 1/2,3/2,5/2, ... (see [57J and [58] for more details). 



(3.7) 



Although a set of useful recurrence relations between the relativistic matrix elements was derived 
by Shabaev [50] , [51] (see also p], [TH], [M], [57], [HI], and the last section of this paper) on the 
basis of a hypervirial theorem, the corresponding relativistic Kramers-Pasternack type relations 
seem to be missing in the available literature. Our equations ( I3.4p - fl3.6l) reveal that they should 
have a vector form. 

Here, one can apply a familiar three-term recurrence relation for the Hahn polynomials 

xhfc to (x, N) = a m h^ (x, N) + f3 m h% « (x, N) + 7m h { «'J ] (x, N) 



with 



a, 



Pr. 



It, 



(m+l)(a + P + m + 1) 
[a + (3 + 2m + 1) (a + (3 + 2m + 2) 



a - (3 + 2N - 2 ((3 2 - a 2 ) (a + (3 + 2N) 

4 + 4 (a + /3 + 2m) (a + /? + 2m + 2) 

(a + m) (13 + m) (a + (3 + N + m) (N - m) 



(a + 13 + 2m) (a + (3 + 2m + 1) 
three times for the corresponding special cases in fl3.4l) - fl3.6l) . Introducing two sets of vectors 

/ -1-2*) \ 



(3.8) 

(3.9) 
(3.10) 
(3.11) 











Bp 


x p - 




\ c P J 





V 



1,(0, 0) 



n-l,-l-2i/ 
(n,l -2i/ 



(3.12) 
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we conclude that the recurrence relations in question should have the following matrix structure: 

Ap = TXp, Xp = T^Ap (3.13) 

with 

Xp+i = UpKp + VpX p _!, (3-14) 
where U p and Wp are known diagonal matrices. Then 



l p+i 



DpA-p -\- EpA-p—i 



(3.15) 
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where 

D p = TUpT' 1 , E p = TVpT- 1 (3.16) 

are similar matrices. 

The elementary linear algebra argument outlined above allows us to obtain a relativistic gener- 
alization of the Kramers-Pasternack recurrence relation (12. 5 j) as follows: 

(2a(3) 2 A p+1 = 4/fe/i^±! A p (3.17) 

p + 2 



pip + 2) (4z/ 2 -(p+l) 2 ) +Ak 



2 

A r 



(p + l)(p + 2) 
V+2^- 1+ (j9+l)(p + 2) ^ 



2p + 3 

(2a^) 2 B p+l = APefx^—- Bp (3.18) 

P ~T~ £ 

-(4z/ 2 -p(p + 2))^5 p _ 1 

p + 2 p + 2 

(2a/3) 2 C p+1 = 4/35/1^1 C p (3.19) 

p + 1 

-Ap-l — ^ — „ -Dp-1 



(p+l)(p + 2) p ^p+2 
p(p + 2) (4z/ 2 - (p + 1) 2 ) -4/i 2 



(p + 1) (p + 2) 

Our relations (I3.17}) - (l3.19j) . together with the symmetry under the reflections, p — > —p — 1 and 
p — > — p — 3, established in [57], give a possibility to compute all the required matrix elements 
in a recursive manner without direct evaluation of the integrals. The corresponding initial values 
(IA.3p - flA.4p can be found in Ref. [57J or elsewhere. 



A different representation of these integrals in terms of Chebyshev polynomials is also available 
[57]: 

2aLi(2a(3) p A p (3.20) 
= a(ji + an) h p ^ 0) (n - 1, -2u) 
+2e (fx 2 - aV) (n - 1, -1 - 2v) 

+a (fu - an) h p ' 0) (n, -2v) , 

2aLi(2ap) p B p (3.21) 
= ae(ii + an) hf' 0) (n - 1, -2z/) 
+2 (/i 2 - a 2 /t 2 ) h^-i } (n - 1, -1 - 2z/) 
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+ae (fj, - an) h p ' 0) (n, -2u) , 

4n(2al3) p C p (3.22) 
= ae(/j, + an) hf' 0) (n - 1, -2v) 
—as (// — an) h p ' ^ (n, —2v) , 
when p > 1. A similar matrix manipulation results in another recurrence relation: 

fo as* A ar (a 2 (P+1) 2 -1) (2p + l) 

2a/? A p+1 = y A p 3.23 

a 2 p (p + 1) (p + 2) 

. 4^V(2p+l) .Jfr+l 
a 2 p (p + l)(p + 2) p p+1 p 
(a 2 (p + l) 2 - 1) (4z, 2 - p 2 ) g (4^ 2 - p 2 ) 



a 2 (p+l)(p + 2) p a 2 (p+l)(p + 2) 

f9^ 2 R 4/feV(2p+l) 

(2a/?) " "a 2 p(p + l)(p + 2)^ (3 ' 24) 



(a 2 (p + l) 2 + g 2 ) (2p+l) _ 2p + l 
a 2 p (p + 1) (p + 2) p p+1 



e(4z/ 2 -p 2 ) (4z, 2 -p 2 )(a 2 (p+l) 2 + g 2 ) 

P _1 9~7 i 1~W T~0~\ -Op— 1) 



a 2 (p+l)(p + 2) p a 2 (p + l)(p + 2) 

(2a/?) 2 C^+i = -^^±1a p + (3s^^ B p (3.25) 

p+1 p+1 

+4/fe/^±i C p - f> 2 - p 2 ) -4- C? H . 
p+1 v p+1 

These recurrence relations are summarized for the benefit of the reader in matrix form in appendix A. 
One should take special values (IA.4I) - (IA.5I) as the initial data. A direct derivation of the relativistic 
three-term vector recurrence relations f)3.17l) - fl3.19l) and (I3.23p - fl3.25l) on the basis of a hypervirial 
theorem needs to be found. 



4. Independent Integrals 

The integrals f)3.1l) - fl3.3p are linearly dependent: 

(2k + e(p+ 1)) A p - (2en + p + 1) B p = 4/jCp (4.1) 

(see, for example, [I], [50], [5T) . and (57j for more details). Thus, eliminating C p , say from (13 . 23ft - 
f)3.24p . we obtain the following three-term vector recurrence relation between the integrals A p and 
B p only: 

(2a/5) 2 A p+l (4.2) 
| ^ 4e/i 2 (a 2 (p + l) 2 - 1) - a 2 p (p + 2) (2k + 5 (p + 1)) 
" (2P + 1} ^ a 2 ^p(p+l)(p + 2) Ap 
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+ (2P + 1] P aV(p+l)(p + 2) Bp 

(a 2 (p + l) 2 -l)(4z/ 2 -p 2 ) e(4z/ 2 -p 2 ) 

P-l 57 i W T~^T -Dp-1 



and 



a 2 (p + l)(p + 2) p a 2 (p + 1) (p + 2) 



(2a/?) 2 (4.3) 

m,, n + fl2 P (g + 2 ) (2k + g (p + 1)) , 
— — (zp -+- ij — - — — — — — — Ap 

a 2 fxp (p + 1) (p + 2) 

4/i 2 (a 2 (p+l) 2 + £ 2 ) +a 2 p(p + 2)(2^ + p+l) 

+ 2p + 1 (3e i f-— — B p 

a 2 fxp (p + 1) (p + 2) 

e(4^ 2 -p 2 ) (4^ 2 -p 2 )(a 2 (p + l) 2 + £ 2 ) 

P _1 2? i\/ T~oN -Dp-1- 



a 2 (p+l)(p + 2) p a 2 (p+l)(p + 2) 

In a similar fashion, from fl3.17l) - fl3.18l) : 



(2a(3) 2 A p+1 = ^sfMj^- A p (4.4) 



2en-(p + 2) (4^ 2 -(p + l) 2 ) 
'"" (p+l)(p + 2) 



and 



2 £K -l + (p + l) (2p + 3) 

— ZK ; r~, ; d„_i 

(p+l)(p + 2) 

9 2n + 3 

(2a/?) 2 5 P+1 = 4/fe/i^— - B p (4.5) 

Matrix forms of these identities and the corresponding initial values are given in appendix A. 

5. Two-Term Recurrence Relations 

The three-term recurrence relations for the relativistic radial integrals, examples of which have 
been found in the previous sections, are obviously not unique. Moreover, if 

A p+1 = D^ + ^Vi, (5-1) 
A p+1 = D^Ap + ^Vi. 

then 

(a p + (3 p )A p+1 = (a p D^ + (3 p D^) A p (5.2) 
+ (a p E^ + /3 p E®) A p _! 
for two arbitrary sequences of scalars a p and (3 p . One special case is of a particular interest. 
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Subtracting (j4.2p - (|4.3p from (I4.4l) - fl4.5p one gets the following matrix equation 

p,(1:)=Q,(1r l ) (5.3) 



with 



Bp I \ B p _i 



, ko 2 (p + 2) (2p + 1) + 2etf 



detP p = -8/3 2 £ ^ — — , (5.4) 

P p(p + l)(p + 2) 2 

detQp = -2e (4z/ J -p 2 ) '-^ — (5.5) 

P K 1 a 2 (p + 2) 2 (p+l) 2 1 ' 

by a computer algebra system. We have omitted the explicit forms of the P and Q matrices and 
their inverses. It should be easily done by the reader. 

The two-term recurrence solutions of the form 



Bp J P I B p _i J 1 B p -i J p I Bp 



where found by Shabaev [5D], [5T] by a different method. In our notations, 

A 4i/ 2 e + 2«(p + 2) + e (p+l)(2« e + p + 2) 

Ap+1 " ~ {P+1 > 4(1 -e 2 ) (p + 2)^ Ap (5 ' 7) 



4/i 2 (p + 2) + (p + 1) (2/ce + p + l) (2/ce + p + 2) 

4(l-e 2 )(p + 2)/?/i p ' 



and 



o r „ 1 n 4^ 2 + 2/ tg (2p + 3)+ g 2 (p+l)(p + 2) 

4/x 2 e (p + 2) + (p + 1) (2k£ + p + 1) (2k + e (p + 2)) 
+ 4(l-£ 2 )(p + 2) / 3/i p 

/3 4/i 2 e(p+l)+p(2/ tg + p)(2/ t + g (p + l)) 
\i (Aw 2 — p z )p 



Ap? (p + 1) + p (2ne + p) (2/te + p + 1) 

^ //I 2 2\ P' 

/i (4z/ z — p z )p 



= ^ + 2..(2 P+ 1) + ^( P+ 1) (5 1Q) 

ji (Av z — p l ) 



4^g + 2k (p + 1) + (2/ce + p + 1) 

^ I A 2 2\ P> 

/i (4^ z — p z ) 



respectively. 



We have shown that these solutions can be derived from the three-term vector recurrence relations 
found in this paper. As a by-product, the factorization of Shabaev's matrices, namely, 

S p = P- 1 Q p , S^ 1 = Q- 1 P p , (5.11) 
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is given. This can be directly verified with the help of a computer algebra system. Then 

det Q p _ (4z/ 2 -p 2 )p 



det S, 



p det P p 4 (a/3) 2 (p+1)' 



(5.12) 



Further details are left to the reader. 
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Appendix A. Matrix Form of Three-Term Recurrence Relations 



The matrix structures of our recurrence relations fl3.17p - fl3.19p and fl3.23p - fl3.25p are given by 



(2apy 




2p + 3 



p + 2 

( p(p+2)(4i^ 2 -(p+l) 2 )+4K 2 
(p+l)(p+2) 



Ap 
Bp 

C 



p+2 



p+2 



2k/x 



(4^-p(p + 2))H±| 



(p+l)(p+2) 



and 



(2a/?) 




2/i 



p+i 

p+2 



p(p+2)(4^ 2 -(p+l) 2 )-V 
(p+l)(p+2) 



p+1 



/ 4e M (a 2 (p+l) 2 -l] 



2p+l 
p+1 



\ 



(Au 2 -p 2 ) 
p+1 

respectively. The initial vectors are 



a 2 p(p+2) 

_ 4e 2 /x 
a 2 p(p+2) 

-1 

/ (a 2 (p+l) 2 ~l) 
a 2 (p+2) 



4e 2 /j 
a 2 p(p+2) 

4efi(a 2 (p+l) 2 +e 2 ) 
a 2 p(p+2) 

£ 



a 2 (p+2) 





a 2 (p+2) 
(a 2 (p+l) 2 + £ 2 ) 

a 2 (p+2) U 

p ) 



-4 \ 



A p 
Bp 



A p -! 

Bp-i 
Cp-i 





I3a 2 

2/11/ • 
1 

e 



( 



) 



(A.l) 




p-i 



(A.2) 



(A.3) 



(A.4) 
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(V 









2^(l-e 2 ) 



\ 



2en-l 
4/3 



(A.5) 



/ 



One should use the initial data (1A.3j) - (lA.4|) in the case of the recurrence relation (1A. 1|) and (1A.4j) - 
fTA~5|) for flA~2|) . 



Our relations fl4.2p - fl4.3l) and fl4.4p - fl4.5l) take the following matrix forms: 

A r . 



(2a(3Y 



X 



(A.6) 



and 



2p + 1 
p + 1 

A p 

B T , 



(2a/?) 1 



I 4eii 2 (a 2 (p+1) 2 -l) -a 2 p(p+2)(2K+e(p+l)) 



— e 



a 2 fip(p+2) 
4ey 2 +a 2 p(p+2) {2n+e (p+ 1) ) 



4eV 2 +a 2 p(p+2)(2eK+p+l) 
a 2 fip(p+2) 

_4/i 2 (a 2 (p+l) 2 +£ 2 )+a 2 p(p+2)(2e K +p+l) 



a 2 fip{p+2) 



(4u 2 - p 2 ) 



a 2 (p + 1) (p + 2) 



a 2 (p+iy-1 

—e 



d 2 iip(p+2) 



a 2 (p+l f + e 2 



A p ^ 
Bp-i 



Ap+i 
Bp+i 



2p + 3 ( A, f 
p + 2 



B r , 



(A.7) 



2£ K -(p+2)(4i/ 2 -(p+l) 2 



n " 'v " ' / Q 2 £K -l+(p+l)(2p+3) \ / A 
P (p+l)(p+2) (p+l)(p+2) ^P- 1 



p+2 



respectively. 



References 

[1] G. S. Adkins, Dirac-Coulomb energy levels and expectation values. Am. J. Phys. 76 (2008) jj=Q, 579-584. 

[2] D. Andrae, Recursive evaluation of expectation values (r k ) for arbitrary states of the relativistic one-electron 

atom, J. Phys. B: Atom. Molec. Phys. 30 (1997) #20, 4435-4451. 
[3] L. Armstrong, Jr., Mean values of powers of the radius for hydrogenic electron orbits, Phys. Rev. A 3 (1971) #5, 

1546-1550. 

[4] R. A. Askey, Continuous Hahn polynomials, J. Phys. A: Math. Gen. 18 (1985), L1017-L1019. 

[5] N. M. Atakishyev and S. K. Suslov, The Hahn and Meixner polynomials of imaginary argument and some of 

their applications, J. Phys. A: Math. Gen. 18 (1985), 1583-1596. 
[6] W. N. Bailey, Generalized Hypergeometric Series, Cambridge University Press, Cambridge, 1935. 
[7] S. Balasubramanian, A simple derivation of the recurrence relation for (r N ), Am. J. Phys. 68 (2000) #10, 

959-960. 

[8] A. I. Baz, Ya. B. Zel'dovich, and A. M. Perelomov, Scattering, Reactions and Decays in Nonrelativistic Quantum 

Mechanics [in Russian], Nauka, Moscow, 1971. 
[9] H. Beker, A simple calculation of (l/r 2 ) for the hydrogen atom and the three-dimensional harmonic oscillator, 

Am. J. Phys. 65 (1997) #11, 1118-1119. 
[10] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms, Springer- Verlag, Berlin, 

1957. 

[11] P. Blanchard, A new recurrence relation for hydrogenic radial matric elements, J. Phys. B: At. Mol. Phys. 7 
(1974) #9, 993-1005. 



RELATIVISTIC KRAMERS-PASTERNACK RECURRENCE RELATIONS 



11 



K. Bockasten, Mean values of powers of the radius for hydrogenic electron orbits, Phys. Rev. A 9 (1974), 1087- 
1089. 

E. Chacon, D. Levi, and M. Moshinsky, Lie algebras in the Schrddinger picture and radial matrix elements, J. 
Math. Phys. 17 (1976) #10, 1919-1929. 

R. Cordero-Soto and S. K. Suslov, Expectation values (r p ) for harmonic oscillator in R n , arXiv: 0908.0032v3 
[math-ph] 6 Aug 2009. 

L. J. Curtis, Classical mnemonic approach for obtaining hydrogenic expectation values r p , Phys. Rev. A 43 
(1991) #1, 568-569. 

L. Davis, A note on the wave functions of the relativistic hydrogen atom, Phys. Rev. 56 (1939), 186-187. 
Yi-B. Ding, On the Schrddinger radial ladder operator, J. Phys. A: Math. Gen. 20 (1987) #18, 6293-6301. 

G. W. F. Drake and R. A. Swainson, Expectation values of r p for arbitrary hydrogenic states, Phys. Rev. A 42 

(1990) #3, 1123-1126. 

J. H. Epstein, S. T. Epstein, Some applications of hypervirial theorems to the calculation of average values, 
Am. J. Phys. 30 (1962) #4, 266-268. 

S. T. Epstein, J. H. Epstein, and B. Kennedy, Recursion relations for Coulomb matrix elements, J. Math. Phys. 
8 (1967) #9, 1747. 

A. Erdelyi, Higher Transcendental Functions, Vols. I— III, A. Erdelyi, ed., McGraw-Hill, 1953. 
A. Gumberidze et ai, Quantum electrodynamics in strong electric fields: The ground state Lamb shift in hydro- 
genlike uranium, Phys. Rev. Lett. 94 (2005), 223001: 1-4. 

A. Gumberidze et ai, Precision tests of QED in strong fields: experiments on hydrogen- and helium-like uranium, 
J. Phys. Conf. Scr. 58 (2007), 87-92. 

J. D. Hey, On the momentum representation of hydrogenic wave functions: Some properties and an application, 
Am. J. Phys. 61 (1993) #1, 28-35. 

J. D. Hey, Further properties of hydrogenic wave functions, Am. J. Phys. 61 (1993) #8, 741-749. 

S. Karlin and J. L. McGregor, The Hahn polynomials, formulas and an application, Scripta Math. 26 (1961), 

33-46. 

S. G. Karshcnboim et al (Eds.), The Hydrogen Atom: Precision Physics of Simple Atomic Systems, Lecture 
Notes in Physics, Vol. 570, Springer- Verlag, Berlin, Heidelberg, 2001. 

S. G. Karshenboim and V. B. Smirnov (Eds.), Precision Physics of Simple Atomic Systems, Lecture Notes in 
Physics, Vol. 627, Springer- Verlag, Berlin, Heidelberg, 2003. 

R. Koekoek and R. F. Swarttouw, The Askey scheme of hypergeometric orthogonal polynomials and its q- 
analogues, Report 94-05, Delft University of Technology, 1994. 

H. T. Koelink, On Jacobi and continuous Hahn polynomials, Proc. Amer. Math. Soc. 124 (1996), 887-898. 

H. A. Kramers, Die Grundlagen der Quantenthcorie/Quantentheoric dcs Electrons und der Strahlung , Akademis- 
che Verlagsgesellschaft, Leipzig, 1938; Engl, transl., Quantum Mechanics, North-Holland, Amsterdam, 1957, 
section 59, equation (6.74). 

L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Nonrelativistic Theory, Pergamon Press, Oxford, 1977. 
R. L. Liboff, Introduction to Quantum Mechanics, third edition, Addison- Wesley, Longman, Reading, MA, 
1998, pp. 475-476. 

H. Marxer, Exact corresponding relationship for the expectation values r~ k for hyprogenlike states, Phys. Rev. 
A 44 (1991) #3, 1543-1548. 

A. Messiah, Quantum Mechanics, Vol. 1, North-Holland, Amsterdam, 1961, p. 431. 

P. J. Mohr, G. Plunicn, and G. Soff, QED corrections in heavy atoms, Phys. Rep. 293 (1998), 227-369. 

R. M. More, Inversion symmetry of hydrogenic expectation values, J. Phys. B: At. Mol. Phys. 16 (1993) #21, 

3859-3861. 

B. Moreno, A. Loperz Pineiro, and R. H. Tipping, Algebraic solutions for the hydrogenic radial Schrddinger 
equation: matrix elements for arbitrary powers of several r-dependent operators, J. Phys. A: Math. Gen. 24 

(1991) #2, 385-390. 

A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov, Classical Orthogonal Polynomials of a Discrete Variable, 
Springer-Verlag, Berlin, New York, 1991. 

A. F. Nikiforov and V. B. Uvarov, Special Functions of Mathematical Physics, Birkhauser, Basel, Boston, 1988. 
H. N. Nunez- Yepez, J. L. Lopez-Bonilla, and A. L. Salas-Brito, Generalized hypervirial and recurrence relations 
for hydrogenic matrix elements, J. Phys. B: At. Mol. Phys. 28 (1995), L525-L529. 



12 SERGEI K. SUSLOV 

[42] P. C. Ojha and D. S. F. Crothers, On a simple relation between hydrogenic radial matrix elements, J. Phys. B: 

At. Mol. Phys. 17 (1984) #24, 4797-4800. 
[43] L. C. Owono Owono, M. G. Kwato Njock, and B. Oumarou, Expectation values of r q between Dirac and 

quasirelavistic wave functions in the quantum- defect approximation, Phys. Rev. A 66 (2002), 052503. 
[44] S. Pasternack, On the mean value of r s for Keplerian systems, Proc. Nat. Acad. Sci. 23 (1937) #2, 91-94; #4, 

250 (erratum). 

[45] S. Pasternack, A generalization of the polynomials F n (x), London, Edinburgh, Dublin Philosophical Magazine 

and J. Science, Ser. 7 28 (1939), 209-226. 
[46] S. Pasternack and R. M. Sternhcimer, An orthogonality property of hydrogenlike radial functions, J. Math. Phys. 

3 (1962) #6, 1280. 

[47] A. M. Perelomov and Ya. B. Zel'dovich, Quantum Mechanics: Selected Topics, World Scientific Publishing Co., 
Inc., River Edje, NJ, 1998. 

[48] W.-C. Qiang and S.-H. Dong, An alternative approach to calculation the mean values r k for hydrogen-like atoms, 

Phys. Scr. 70 (2004), 276-279. 
[49] M. L. Sanchez, B. Moreno, and A. Loperz Pineiro, Matrix-element calculations for hydrogenlike atoms, 

Phys. Rev. A 46 (1992) #11, 6908-6913. 
[50] V. M. Shabaev, Recurrence formulas and some exact relations for radial integrals with Dirac and Schrddinger 

wave functions, Vestnik Leningradskogo Universiteta, Seria 4: Fizika, Khimiya (1984) number 1, 15-19 [in 

Russian] . 

[51] V. M. Shabaev, Generalizations of the virial relations for the Dirac equation in a central field and their appli- 
cations to the Coulomb field, J. Phys. B: At. Mol. Opt. Phys. 24 (1991), 4479-4488. 

[52] V. M. Shabaev, Relativistic recoil corrections to the atomic energy levels, in: The Hydrogen Atom: Precision 
Physics of Simple Atomic Systems, Lecture Notes in Physics, Vol. 570, Springer- Verlag, Berlin, Heidelberg, 
2001, pp. 714-726. 

[53] V. M. Shabaev, Two-time Green's function method in quantum electrodynamics of high-Z few-electron atoms, 

Phys. Rep. 356 (2002), 119-228. 
[54] V. M. Shabaev, Virial relations for the Dirac equation and their applications to calculations of hydrogen-like 

atoms, in: Precision Physics of Simple Atomic Systems, Lecture Notes in Physics, Vol. 627, Springer- Verlag, 

Berlin, Heidelberg, 2003, pp. 97-113. 
[55] V. M. Shabaev, Quantum electrodynamics of heavy ions and atoms: current status and prospects , Phys.- Uspekhi 

178 (2008) #11, 1220-1225 [in Russian]; Phys.- Uspekhi 51 (2008) #11, 1175-1180. 
[56] J. Shertzer, Evaluation of matrix elements (n, I \r@ | n, l ,s f for arbitrary [3, Phys. Rev. A 44 (1991) #5, 2832-2835. 
[57] S. K. Suslov, Expectation values in relativistic Coulomb problems, J. Phys. B: Atom. 

Mol. Opt. Phys. 42 (2009), 185003 (8pp); pub lished on line 9 September 2009, URL: 

|http://www.iop.org/EJ/abstract/0953-4075/42/18/185003/[ 
[58] S. K. Suslov and B. Trey, The Hahn polynomials in the nonrelativistic and relativistic Coulomb prob- 

lems, J. Math. Phys. 49 (2008) #1 , 012104 (51pp); published on line 22 January 2008, URL: 

|http://link.aip.org/link/?JMP/49/012104[ 
[59] R. A. Swainson and G. W. F. Drake, An alternative proof of some relations between hydrogenic matrix elements, 

J. Phys. B: At. Mol. Opt. Phys. 23 (1990) #7, 1079-1084. 
[60] P. L. Tchebychef, Sur V interpolation par la methode des moindres carres, Memoires de l'Academie Imperiale des 

sciences de St.-Petersbourg, VIP serie, 1 (1859) #15, 1-24; reprinted in: Qeuvres, volume 1, Chelsea Publishing 

Company, New York, 1962, pp. 473-498. 
[61] P. L. Tchebychef, Sur V interpolation, (1864); reprinted in: Qeuvres, volume 1, Chelsea Publishing Company, 

New York, 1962, pp. 542-560. 
[62] P. L. Tchebychef, Sur V interpolation des valeurs equidistantes , (1875); reprinted in: Qeuvres, volume 2, Chelsea 

Publishing Company, New York, 1962, pp. 219-242. 
[63] J. H. Van Vleck, A new method of calculating the mean value of l/r s for Keplerian systems, Proc. Roy. 

Soc. London A 143 (1934), 679-681. 
[64] E. R. Vrscay and H. Hamidian, Rayleigh-Schr 6 dinger perturbation theory at large order for radial relativistic 

Hamiltonians using hypervirial Hellmann-Feynman theorems, Phys. Lett. A 130 (1988) #3, 141-146. 
[65] I. Waller, Der Starkeffekt zweiter Ordnung bei Wasserstoff und die Rydbergkorrektion der Spektra von He und 

Li+, Zeitschrift fiir Physik 38 (1926), 635-646. 



RELATIVISTIC KRAMERS-PASTERNACK RECURRENCE RELATIONS 



13 



School of Mathematical and Statistical Sciences and Mathematical, Computational and Mod- 
eling Sciences Center, Arizona State University, Tempe, AZ 85287-1804, U.S.A. 

E-mail address: sks@asu.edu 

URL: http : //hahn . la . asu . edu/~suslov/ index . html 



